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It is well established that lattice artifacts can be suppressed substantially by the use of SU(3)- 
projected smeared links in the fermion action. An example is the Highly Improved Staggered 
Quark action where the ASQ-like effective links are constructed from reunitarized Fat7 links. A 
general procedure is presented for computing the derivative of the fermion action with respect 
to the base links (fermion force) — a key component in dynamical simulations using molecular 
dynamics evolution. The method is iterative and can be applied to actions with arbitrary levels 
of smearing and reunitarization. The cost of calculating the fermion force is determined for the 
ASQ action and the HISQ action. Test results show that calculating the HISQ force is about two 
times more expensive than the ASQ force. 
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1. Introduction 

The staggered quark formalism provides a fast method to simulate light quarks at small quark 
masses. Staggered fermions, however, have the property that each lattice quark field describes 
four identical quark flavors ("tastes") rather than one, and coupling to gauge fields leads to taste- 
changing errors, which are 0(a 2 ) effects [p. Taste-changing interactions can be suppressed sub- 
stantially by the use of fat links [jlj, ^, ^, |J ; accurate simulation results [6], [7[] were obtained 
using the ASQTAD action, where the base links are replaced by Fat7 effective links (see Fig. 1). 
The success of the ASQTAD action suggests that further improvement may be possible with more 
smearing. Perturbative calculations and non-perturbative results of Refs. H [J show that signifi- 
cant improvement can be obtained by first projecting the fat links back to SU(3) matrices before 
being used in the next fattening process. The result is the Highly Improved Staggered Quark action 
(HISQ) where the ASQ-like effective links are constructed from reunitarized Fat7 links (another 



widely used smearing technique is HYP-smearing [JlOp , which is based on SU(3)-projected fat links 
also) 

U eff (U R ,U RJf ), U R {U F ,U n ) } U F {U,U^). (1.1) 

Here U is the base link, U F is the Fat7 link, U R is the reunitarized fat-link and U e f f is the ASQ- 
like effective link (see Fig. 1). The HISQ action has no tree-level order a 2 errors like the ASQTAD 
action, and has taste-changing interactions that are 3-4 times smaller than the ASQTAD action. In 
addition, the action has no tree-level order (am) 4 errors to leading order in the quark's velocity v/c, 
it therefore provides an accurate discretization of the charm quark on the lattice. For example this 
action has been used recently to obtain high precision results for D meson and decay constants for 



71, K, D and D s [g|llj|. 

So far all unquenched simulations using staggered quarks as light quarks were done with the 
ASQTAD action. Given the nice features of the HISQ action, it is desirable to use it for the sea 
quarks also. Dynamical simulations with HISQ, however, are complicated by the extra level of 
fattening and SU(3)-projection. In this paper we present a general procedure for computing the 
derivative of the fermion action with respect to the base links (fermion force) — a key component 
in dynamical simulations using molecular dynamics evolution. This method is iterative and can be 
applied to actions with arbitrary levels of smearing or SU(3)-projection. We compare the efficiency 
of the algorithm for the ASQTAD and HISQ actions on small lattices. 



Our treatment of unitarized links have been influenced by Kamleh et al. [12]. Other approaches 
have been discussed in Refs. [T3|, 14, 15]. 



2. The Fermion Force 



The staggered quark action is (we follow the notations of [16]) 

-n f /4 



[M f [U]M[U}]~ 



(2-1) 



where 

M X!y [U] =2m8 Xd , + D x , y [U] 



2md x , y + £ i lxfl [UZV 3 X , „ - UZ%8 W J . (2.2) 
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Figure 1: Paths used in the ASQTAD and HISQ actions. Path coefficients for ASQTAD can be found 
in Ref. [[jj. The HISQ effective links are constructed by first applying a Fat7 fattening to the base 
links (U -> U F with coefficients l-link:l/8, 3-staple:l/16, 5-staple: 1/64, 7-staple: 1/384), then a SU(3)- 
projection (U F — > t/ R ), and finally an "ASQ" smearing (U R — ► with coefficients 1-link: 1 +e/8, 3- 
staple:l/16, 5-staple: 1/64, 7-staple: 1/384, Lepage: — 1/8, Naik: — (1 +e)/24, where the parameter e is in- 
troduced to remove (am) 4 errors [pfl). 



The pseudo-fermion field <J> is defined on even lattice sites only to avoid a doubling of flavors from 
using M^M instead of M in the action. This procedure is valid since M^M has no matrix element 
connecting even and odd lattice sites. 

A key component in dynamical simulations using molecular dynamics evolution is the com- 
putation of fermion force — derivative of the fermion action with respect to the base links 

dS f _ d 



fx,[i 



[M f [U]M[U]] 



-nf/4 



4> 



(2.3) 



The derivative can be computed straightforwardly if tif is a multiple of 4; for other numbers of 
fermion flavors the 4th-root of M^M can be approximated by a rational expansion (the RHMC 



algorithm [17, mn 



[M t M]-"// 4 «ob + £- 



a/ 



,'MtM + jS/ 

where a/ and /3/ are constants. The derivative becomes 



(2.4) 



_dSf_ 
dU x , 



i-i 



(M f [U]M[U]) 



dDi[U] 



dU x , 



Y l ) + (Y' 



dD[U] 



dU x , 



X 1 



(2.5) 



with \X l ) = [M t M + j8/] _1 |<J>) and \Y l ) = D\X l ). Note that X 1 and Y l are defined on even and odd 
sites respectively. Taking the derivatives of D, D T with respect to U'ff, U ei ^ and writing out the 
matrix indices we have 



[f*s)ab 

c(0) 



dS f 



V( IVti I d\Uyy\mn f (0) 

1 / ' j,v 3777 i Jy,v 



+ 



j',v Jmn 



where f^'J is the vector outer product of the field variables at y and y + V 

' £ a/ [Fy+ V ]„ [Xy'*] m for even y 

^ a l\ X l+vWy*]m for odd y ' 



\i\ab 



>(0)t 
Jy,v 



(2.6) 



>(0) 
V-y.v 



(2.7) 
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and the sum on y, v extends over the effective links that contain the base link U Xtfl , and (— l) y = 1 
for even y's and (— l)- y = —1 for odd sites. Note that the force is equal to /(°) when there is no 
smearing, i.e., U e ff = U. Furthermore all complications associated with the rational expansion 
have been absorbed into /(°) and the derivatives dU e f f jdJJ , dU e f^ jdJJ are calculated only once. 



3. Smearing and SU(3)-Projection 

For actions with multiple levels of smearing or SU(3)-projection, such as the HISQ action, we 
use the chain rule to compute dU e ^ /dll, dU e f^/dU 



[fx,fi. 



ah 



I ( 

y,v;z,p 



f ( d[U; f /] mn d[U« p } pq 

\ \d[u* p ] pq d[u x ^] ab 
+ ( u eff -» u efft 



+ 



d[uS f ] mn d[u«X q 



x,\i\ab 



f 



(0) 



f (0)t 
Jy,v 



>(i) 

Jz,p 



PI d[U x ^] a b 



f (l)t 

Jz,p 



(3.1) 



PI . 



where the sum on z, p extends over the reunitarized links U R that contain the base link U x ^, and 



f (i) 



PI 



E(-1F%; 



'J IKV. |>(0)1 , dWtP] 

Jy,v 



z,p)pq 



+ ■ 



z,pip<i 



>(0)t 
Jy,v 



(3.2) 



The expression for f x4i still contains the composite derivative dU R /dU. Repeat this step until all 
the derivatives are explicit 



[fx.jl 



ah 




(3.3) 



Note that / (2) is local since U R p is a function of U F p, U^p only. Therefore to construct the fermion 
force one starts with f@' , computes /W for each smearing and SU(3)-projection working in towards 
the first level of smearing. This procedure is very general and can be applied to actions with 
arbitrary levels of smearing and reunitarization. 

Two types of derivatives are involved: derivatives of smeared links (dU e f f /dU R and dU F /dU) 
and derivatives of reunitarized links (dU R /dU F ). Implementing the smeared links and their deriva- 
tives is relatively straightforward. For SU(3)-projected smeared links we use polar decomposition 
since it is differentiable and therefore suitable for dynamical simulations (another popular choice 



is the stout link introduced by Morningstar and Peardon [13]) 



u 



jr—x, U R = U F [U^U 



Fl-1/2 



(dett/*) 1 / 



(3.4) 
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There are several ways to compute the derivative dU R /dU F . Refs. [|13|, [15|] utilize the Cayley- 
Hamilton Theorem; in Ref. [0] the derivative is obtained by solving a matrix equation. Here we 

—1/2 

follow Ref. [12] by adopting a rational expansion for [£/ F ^£/ F ] 



(3.5) 



where c\ and di are constants. A nice feature of the rational approximation is that the derivatives of 
U R are calculated easily 



d[U R ' 



d[U> 



->pr 



ci 



5> 



u 1 



j<U Ff U F + d, 
1 



sq 



U pJ fU F + di 



pr 



1 



U F W F + d t 



(3.6) 



sq 



and including the determinant 



dp 1 



i pi 



dp 1 



(detU 



- 3 tr[U 



, dU R \ d[U R ~ 



d[U F ] rs 



pi 



dp 1 



(3.7) 



where the trace is taken with respect to the indices on U R . Finally the derivative can also be 
computed numerically 



dp%q _ p R ] Pq {p F ] rs + h)-p R ] pq (p F ] rs -h) 

dp F ] rs 2h 



(3.8) 



with £/ Ft being held fixed. The matrices U R {U F ±h) can be obtained by diagonalizing U F U F ^. 
This method, however, is relatively inefficient compared to the rational approximation approach 
(typically by a factor of five). On the other hand, it provides a useful check on our approximate 
calculation. 



4. Benchmark 

We test the performance of the HISQ force on small lattices with a scalar code. We first 
compare the average plaquettes obtained using the RHMC algorithm and the R-algorithm. A high 
acceptance rate provides an excellent check of the code. Results are shown in Fig. 2. 

Fig. 3a compares the cost of different components for the ASQTAD action and the HISQ action 
at different lattice volumes. Results show that matrix inversions (with no optimization) dominate 
for ASQTAD while the computation of fermion force and inversions are comparable for HISQ. It 
is interesting to note that the cost of calculating f&\ derivatives of reunitarized links, is relatively 
small compared to /W and /( 3 ), derivatives of smeared links. This implies that the HISQ force is 
only about twice as expensive as the ASQ force. This point is further emphasized in Fig. 3b where 
the ratio of the cost of fermion force to the cost of gauge force is plotted for the two actions. Finally, 
the insensitivity of our results to lattice volume suggests that Fig. 3 is valid for larger lattices also. 
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Figure 2: Average plaquettes obtained using the RHMC algorithm and the R-algorithm for the HISQ action 
(with e = 1). Simulation parameters are j3 p / = 8.0, rif = 2, am = 0.25, V = 4 4 and n„ K / x dt ^ 0(1) where 
dt = 0.02,0.04,0.06 is the molecular dynamics step size and n mc i = 50,25, 16 is the number of steps. The 
1-loop Symanzik-improved gluon action is used with tadpole factor uq = 0.8897. The RHMC acceptance 
rates are shown in red. 




■ f(3) gauge icomp. of ot h e rs L 

force eff. links 



Figure 3: a) Cost of different components: /(°) is the force when there is no smearing, /«and/( 3 ) are the 
derivatives of smeared links and /' 2 ' is the derivative of reunitarized links. Simulation parameters are the 
same as those in Fig. 2 and the RHMC algorithm is used, b) Ratio of the cost of fermion force to the cost of 
gauge force for ASQTAD and HISQ. 
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5. Conclusion 

A general procedure is presented for computing the fermion force for actions constructed with 
SU(3)-projected smeared links. Application to the HISQ action has been discussed and tests have 
been done on small lattices. Our results show that the HISQ force is only two times more expensive 
than the ASQ force, with most of the additional cost attributed to the extra level of Fat7-smearing. 
Given that it is relatively inexpensive to compute the fermion force one should seriously consider 
using SU(3)-projected smeared links in future simulations. 
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